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ABSTRACT 

Rapid magnetic upflows in the quiet-Sun photosphere were recently uncovered from both SUNRlSE/IMaX and 
Hinode/SOT observations. Here, we study magnetic upflow events (MUEs) from high-quality, high- (spatial, 
temporal, and spectral) resolution, and full Stokes observations in four photospheric magnetically sensitive 
Fe I lines centered at 5250.21, 6173.34, 6301.51, and 6302.50 A acquired with the Swedish Solar Telescope 
(SST)/CRISP. We detect MUEs by subtracting in-line Stokes V signals from those in the far blue wing whose 
signal-to-noise ratio (S/N) > 7. We find a larger number of MUEs at any given time (2.0 x 10 -2 arcsec -2 ), 
larger by one to two orders of magnitude, than previously reported. The MUEs appear to fall into four classes 
presenting different shapes of Stokes V profiles with (I) asymmetric double lobes, (II) single lobes, (III) double¬ 
humped (two same-polarity lobes), and (IV) three lobes (an extra blueshifted bump in addition to double lobes), 
of which less than half are single-lobed. We also find that MUEs are almost equally distributed in network and 
internetwork areas and they appear in the interior or at the edge of granules in both regions. Distributions of 
physical properties, except that of horizontal velocity, of the MUEs (namely, Stokes V signal, size, line-of-sight 
velocity, and lifetime) are almost identical for the different spectral lines in our data. A bisector analysis of 
our spectrally resolved observations shows that these events host modest upflows and do not show a direct 
indication of the presence of supersonic upflows reported earlier. Our findings reveal that the number, types 
(classes), and properties determined for MUEs can strongly depend on the detection techniques used and the 
properties of the employed data, namely, S/Ns, resolutions, and wavelengths. 
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1. INTRODUCTION 

The solar photosphere is a highly dynamic environment 
where motions with a wide range of spatial and temporal 
scales are observed (Spruit et al. 1990). Both magnetized and 
non-magnetized flows are found in the photosphere. The mag¬ 
netized flows are related to or driven by the magnetic fields 
whose manifestation in the solar photosphere are on all spa¬ 
tial scales (Stenflo 1989; Solanki 1993, 2001; de Wijn et al. 
2009). Among those, observations of small-scale magnetic 
fields and their interactions with convective motions at high 
spatial and temporal resolutions have uncovered a number of 
interesting phenomena in the solar photosphere (e.g., Bellot 
Rubio et al. 2001; Shimizu et al. 2007; Fischer et al. 2009; 
Martinez Gonzalez & Bellot Rubio 2009; Zhang et al. 2009; 
Danilovic et al. 2010; Steiner et al. 2010; Straus et al. 2010). 
Shimizu et al.’s (2008) discovery of high-speed (supersonic; 
Fd 8 -10 kms -1 ) downflows associated with the formation of 
small concentrated magnetic patches is one example of those 
small-scale magnetized flows to particularly mention here. 

The magnetic fields in the solar atmosphere are inferred 
from measurements of polarization states of light through ob¬ 
servations of magnetically sensitive lines (Wittmann 1974; 
Auer et al. 1977; see also Borrero & Ichimoto 2011 and refer¬ 
ences therein). The polarization states are represented by four 
Stokes parameters denoted by I referring to the total intensity, 
Q and U describing linear polarizations, and V characterizing 
circular polarization (Unno 1956; Rachkovsky 1962; Stenflo 
1971; Landi DeglTnnocenti & Landolfi 2004). 

Recently, Borrero et al. (2010) claimed discovery of su¬ 
personic magnetic upflows from strong Stokes V signals they 
found in the continuum position of the IMaX Fe I 5250.2 A 
passband (Martmez-Pillet et al. 2011) aboard the Sunrise 
balloon-borne solar observatory (Solanki et al. 2010). Those 


large velocities were, however, inferred based on a single con¬ 
tinuum position, from which the shape of Stokes V profiles 
could not be sampled. The existence of these magnetic up¬ 
flows was later confirmed by Martinez Pillet et al. (201 1) who 
referred to them as photospheric jets. They made use of data 
from Hinode/SP (Tsuneta et al. 2008), which have lower spa¬ 
tial resolution, but a higher spectral sampling compared to the 
data from Sunrise. Using the IMaX data, Rubio da Costa 
et al. (2015) showed the number density of the events is in¬ 
dependent of the heliocentric angle. They found a number 
density of ^ 7.0 x 10 -4 arcsec -2 , which is smaller than the 
values found by Borrero et al. (2010) and Martinez Pillet et al. 
(2011) (« 2.0 x 10 -3 arcsec -2 ), but comparable to the number 
density reported by Borrero et al. (2013). 

These events were explained in terms of magnetic recon¬ 
nection, where freshly emerged flux tubes within granules in¬ 
teract with pre-existing fields in intergranular areas (Borrero 
et al. 2013; Quintero Noda et al. 2013). Quintero Noda et al. 
(2014) found that the extremely blueshifted events, observed 
by Hinode/SP, are often accompanied by redshifted Stokes V 
signals. Together, they represent footpoints of magnetic U- 
loops with flows along them. Quintero Noda et al. (2014) ar¬ 
gued that the siphon flow mechanism (Montesinos & Thomas 
1993) was the only explanation for those flow motions along 
the arched magnetic flux tubes. 

On the other hand, in the MHD simulations of Danilovic 
et al. (2015), the rapid magnetic flows could solely be pro¬ 
duced by flux emergence (e.g., Caligari et al. 1995; Lites et al. 
1996; Orozco Suarez et al. 2008). While magnetic reconnec¬ 
tion could also take place in these simulations, it was not, nec¬ 
essarily, the cause of such events. 

In the above studies, only single-lobed Stokes V profiles 
were reported at the spatial locations of the rapid magnetic 
upflows, albeit Martinez Pillet et al. (2011) also noted ob- 
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servations of highly asymmetric Stokes V profiles with mul¬ 
tiple lobes (shifted with respect to the rest wavelength po¬ 
sition). However, they did not analyze the profiles. We 
should note that few cases of supersonic magnetic upflows 
with double-peaked structures in the blue lobe of Stokes V 
profiles (i.e., three-lobed Stokes V profiles, where the blue 
lobes are double-peaked) were reported previously by Socas- 
Navarro & Manso Sainz (2005). 

In this paper, we detect and study a wide range of magnetic 
Doppler shifted phenomena and in particular magnetic upflow 
events (MUEs) from observations with not only high spatial 
and temporal resolutions but also with highly sampled profiles 
(from the Swedish Solar Telescope (SST)/CRISP; Section 2). 
We also conduct our investigations on a time-series of images 
sampled in four different photospheric lines. In addition, we 
use a different method compared to that employed in previous 
studies with which a larger number of events (larger by one or 
two orders of magnitudes) are identified (Section 3). Distri¬ 
butions of the physical and dynamical properties of the events 
are studied and their spatial distribution over network and in¬ 
ternetwork areas is inspected (Section 4). We find four differ¬ 
ent classes of Stokes V profiles corresponding to the MUEs. 
Finally, we discuss our results, compare them with those in 
the literature, and draw our conclusions in Section 5. 

2. DATA 

In the present work, we use several high-quality, high- 
resolution data sets recorded at the quiet-Sun disk center 
by the CRisp Imaging SpectroPolarimeter (CRISP; Scharmer 
et al. 2008) mounted at the SST (Scharmer et al. 2003). The 
combination of high overall transmission (with a minimum of 
optical surfaces) of the CRISP dual Fabry-Perot Interferome¬ 
ter (FPI) and a high-cadence camera system provide a fast line 
sampling (with a typical effective sampling speed of 0.25 s 
per line position for non-polarimetric observations and 1-1.5 s 
per line position for the polarimetric data analyzed here) for a 
wide wavelength range of 5100-8600 A with numerous pho¬ 
tospheric and chromospheric spectral diagnostics (Scharmer 
2006). In addition, the combination of high and low spectral 
resolution etalons at the CRISP FPI results in a high spec¬ 
tral resolution at a selected wavelength. Polarization mea¬ 
surements are performed by two liquid crystal modulators in 
combination with a polarizing beam splitter which facilitates 
a large reduction of seeing-induced polarization crosstalk by 
simultaneously exposing two cameras. 

The employed full Stokes (/, Q , U, and V) time series 
of images were acquired in four photospheric, magnetically 
sensitive, Fe I lines. Table 1 summarizes characteristics of 
these data sets along with the line parameters. To estimate the 
heights of formation of the four Fe I lines, we computed the 
line depression contribution functions using the RH radiative 
transfer code of Uitenbroek (2001) in FAFC, FAFF, and FAFP 
model atmospheres (Fontenla et al. 1993, 2006). These atmo¬ 
spheric models represent an averaged quiet-Sun area, bright 
regions of the quiet-Sun, and a typical plage region, respec¬ 
tively. We determined the contribution of the atmospheric 
heights to the radiation passing through the line core of our 
SST/CRISP Fe I lines (where the spectra were convolved with 
the transmission profiles of the filters) in non-FTE conditions. 
The formation heights given in Table 1 are the average values 
(weighted by the corresponding contribution functions) from 
the calculations in the FAFC mode atmosphere. The average 
of the formation heights is slightly (less than 10 km) smaller 


Table 1: Overview of Fines and Data Sets from SST/CRISP 


Quantity 

Fe I 5250 A 

Fe I 6173 A 

Fe I 6301/6302 A 

Line center 

5250.21 A 

6173.34 A 

6301.51/6302.50 A 

Lande factori 

3.0 

2.5 

1.67/2.5 

log (gifik) h 

-4.938 

-2.880 

-0.718/-0.973 

Spectral resolution 0 

26.3 mA 

50.7 mA 

54.8 mA 

Formation height d 

214 km 

174 km 

178/125 km 

Date of observations 

2014 Jun 23 

2012 Jun 4 

2011 Aug 6 

Duration 

13 minutes 

41 minutes 

47 minutes 

Cadence e 

20 s 

14 s 

28 s 

Acquisition f 

20 s 

14 s 

18 s g 

No. of positions 11 

21 

10 

15/15 

Sampling step 

25 mA 

35 mA 1 

44 mA 

No. of exposures! 

8 

11 

5 

M k 

1 

1 

1 


Notes. 

a Magnetic sensitivity of the lines. From Socas-Navarro (2004). 
b A proxy for line strength between the lower (/) and upper (k) energy lev¬ 
els where /# is the oscillator strength and g is the statistical weight. From 
Fabbian et al. (2012) and Kramida et al. (2014). 

c FWHM of the CRISP transmission profile. From de la Cruz Rodriguez 
et al. (2015). 

d Average heights of formation of the line cores for FALC model atmo¬ 
sphere, after convolution with the CRISP transmission profiles. 
e Time step in the time series of images. 
f Total time for acquisition of a full spectral scan. 

g This data set also includes observations of another line scan (Fe I 5576 A). 
h Number of wavelength positions. 

1 The width of spectral sampling step between the second and eighth wave¬ 
length positions of this line. It is 70 mA for the first and the ninth and 140 mA 
for the last position. 

j Number of exposures per polarization state. 
k The cosine of the observing angle. 

in the other two atmospheric models compared to the FAFC, 
for all the lines except for the 5250 line. The latter is smaller 
by « 50 km in FAFP compared to that obtained in the FAFC 
model atmosphere. 

The images have a scale of « 0.06 arcsec pixel -1 cover¬ 
ing an average field of view (FOV) of 52 x 52 arcsec 2 after 
restoration. The 6301/2 A data set also includes simultaneous 
observations in the Ca II H passband with a bandwidth of 1.0 

A. 

All recorded images were processed using the CRISPRED 
reduction pipeline (de la Cruz Rodriguez et al. 2015) 
which also restores images from wavefront aberrations by 
means of Multi-Object Multi-Frame Blind Deconvolution 
(MOMFBD; van Noort et al. 2005). The pipeline also 
includes calibrations for field-/time-dependent instrumental 
polarization, inter-camera and temporal misalignments, as 
well as image field rotation, and corrections for reflectivity 
and cavity errors (induced in the CRISP FPI passband pro¬ 
files). We use multiple exposures per polarization state in the 
MOMFBD restoration which increases the signal-to-noise ra¬ 
tio (S/N) (van Noort & Rouppe van der Voort 2008). Ulti¬ 
mately, the number of exposures per state is a trade-off be¬ 
tween the necessity of high S/N dense spectral sampling, and 
the requirement to complete the line scan within the typi¬ 
cal evolution timescale in the photosphere. The adaptive op¬ 
tics system at the SST (Scharmer et al. 2003) as well as the 
MOMFBD post-processing approach result in images with 
high spatial resolution close to the diffraction limit of the 1 
m telescope (\/D = 0'.'l at 5250 A). 

Figure 1 illustrates mean Stokes I profiles of all the pass- 
bands (averaged over the entire FOV) along with their sam- 
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reddest points for the number of corresponding marked circles 
(n\, ri 2 , n~, and n + ). 

The Stokes parameters have a lcr noise level of « 1.6 x 
1CT 3 I c per wavelength, after post-processing. The noise level 
was computed as the standard deviation at a far continuum 
position of the Stokes profiles where almost no signal is ap¬ 
parent. We note that the la noise level of V“ ing , V^ ing and Vn n& 

images are reduced down to (6.0- 11.0) x 10 -4 I c depending 
on the number of wavelength positions used to form the aver¬ 
aged images (see Figure 1). 

For the sake of simplicity and of consistency, we will show 
all examples (images and profiles) in the following for those 
recorded in the Fe I 6301/6302 lines. The same analyses are 
employed for all four spectral lines (unless otherwise stated). 
The number of detected events from the different passbands 
(and some of their physical properties) are, however, different 
(see Section 4). 


Wavelength (A) Wavelength (A) 


3. EVENT DETECTION 


Figure 1.: Averaged Stokes I profiles of the four photospheric 
Fe I lines centered at 5250.21 A (a), 6173.34 A (b), 6301.51 A 
(c), and 6302.50 A (d) whose Stokes V profiles are analyzed 
in this paper. The circles indicate the sampled wavelength 
positions. The green, blue, and red circles (with cross, minus, 
and plus signs) mark those wavelength positions from which 
the averaged in-line, far-blue-wing, and far-red-wing Stokes 
V images used for event detection are constructed. 


pled wavelength positions marked with circles. For the pur¬ 
pose of event detection, we combine several line positions 
to represent different parts of the profile at higher signal-to- 
noise. Different colors/shapes of the circles indicate the se¬ 
lected wavelength positions from which averaged Stoke V sig¬ 
nals are combined in the far blue wing (V“ ing ; blue circles), far 

red wing (V^ ing ; red circles), and in-line positions close to the 
line-core (Vjj ne ; green circles). To avoid cancellation, the signs 
of the wavelength points on the red side of the line core are 
reversed prior to forming the in-line images. These averaged 
Stokes V images have a larger S/N compared to the individual 
wavelength points (the noise levels are reduced by the square 
root of the number of wavelength positions with which the 
averaged images are formed). Hence, their absolute values 
are being used for detection of the magnetic Doppler-shifted 
events (see Section 3). Thus, three sets of images, from each 
Stokes V profile, are formed as 


Inline I = 


1 


n 1 +M 2 


«1 «2 

y>(A;)-y>(A;) 


( 1 ) 




wing | 


1 

nr 


y> A ‘-) 


( 2 ) 


We aim to study Doppler-shifted (in particular blueshifted) 
events associated with the magnetic fields. To this end, we 
want to detect Stokes V signals that are largely shifted toward 
the blue or toward the red wavelengths, i.e., where the Stokes 
V signal in a pixel is significantly larger in the far wings com¬ 
pared to the in-line signal of the same pixel. In the following, 
we refer to these features as magnetic upflow events (MUEs) 
and magnetic downflow events (MDEs). We will, however, 
exclude those events that are smaller than the spatial resolu¬ 
tion of our observations (i.e., ^10 pixels in area). In addition, 
in order to secure detections of true signals, we exclude pix¬ 
els whose far-wings Stokes V signals are smaller than their 
corresponding la noise levels (i.e., 0.6% I c on average; in the 
range of 0.4% -0.8% I c for the different sampled lines). 

For the event detection, we propose an approach that is 
based on subtraction of |Vii ne | from |I/T ng |. To compare our 
method with that employed in the earlier works by, e.g., Bor- 
rero et al. 2013, we also apply a technique that is defined as 
the quotient of the far wings and in-line signals. 

3.1. Difference Method 

Positive values of the difference between absolute values 
of far-wings Stokes V signals (far blue wing, |V“ ing |, and/or 

far red wing ,|V^ ing |) and those of in-line ones (|Vime|) rep¬ 
resent Stoke V signals that are largely shifted toward blue or 
red wavelengths corresponding to magnetic upflows and mag¬ 
netic downflows, respectively. This builds the base for our 
event detection approach. We should, however, make sure that 
the real signal is properly distinguished from the noise prior 
to the subtraction. Thus, we set the la noise level thresh¬ 
old on the |V^ ng | images. The definitions of magnetic up- 
flows/downflows are then denoted as 


ft 


+ | 
wing | 


1 

n + 


n + 


rlef i i 

MUE := |V- ing | sig -|V,ine| > 0, 

(4) 

y>( a, + ) 

i + 

, (3) 

MDE:= |V+ ng | sig -|Vline|>0, 

(5) 


where the indices i, T, and i + represent the scanning wave¬ 
length positions indicated in Figure 1 by the green circles 
with cross signs, blue circles with minus signs, and red cir¬ 
cles with plus signs, respectively, running from the bluest to 
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As noted earlier, we only keep those events whose areas 
are larger than 10 pixels. Finally, we do a visual inspection 
of all individual frames to exclude possible false (or unclear) 
detected events. The latter could be a result of, e.g., im¬ 
age artifacts (fringes) coinciding at the location of the can¬ 
didate MUEs/MDEs. We note that the absolute values of 
MUEs/MDEs signals are only used for detection purposes. 
Both events have been found to be of either polarity. 

3.2. Quotient Method 

For comparison reasons, we also detect MUEs (from the 
same data sets as in Section 3.1) by employing a similar ap¬ 
proach as described by Borrero et al. (2013). In this method, 
the absolute value of the quotient of far-blue-wing and in-line 
Stokes V signals is determined and the MUEs are defined as 
l^wing/^linel ^ 1- 

For a direct comparison with our difference approach, we 
implement the quotient method after masking |V“ ing | with 
2 cr no i S e of the employed data sets, and similarly exclude MUEs 
smaller than 10 pixels. Also, we redo the quotient method 
when |V“ ing /Viine| > 4 and without the earlier criteria, i.e., 
same as employed by Borrero et al. (2013). See Section 5 
for comparisons and discussions. 

3.3. A Case Study 

Illustrated in Figure 2 are examples of the far-blue-wing 
(top left), far-red-wing (bottom left), and in-line (top right) 
Stokes V images along with an image of Stokes I continuum 
(bottom right) recorded in a scan of the 6302 A line (see Sec¬ 
tion 3.1 for the definitions). For better visibility and for fur¬ 
ther visual comparisons, all the Stokes V images are saturated 
to the same scale. The actual signals of the V“ ing , V+ ing , and 
Vii ne images in Figure 2 are in the range of (-2.4%,3.3%), (- 
3.6%,3.0%), and (-7.3%,8.5%), relative to the Stokes I con¬ 
tinuum (7 C ), respectively. The blue circles in Figure 2 include 
the detected MUEs from our difference method. We remind 
the reader that the events may have either polarities and ab¬ 
solute values of their signals are only used to facilitate their 
detection. One of the MUEs is magnified within its circle 
(all open circles have an inner diameter of « 1 arcsec). A vi¬ 
sual comparison of the events, indicated by the blue circles, 
shows a significantly larger Stokes V signal in the far-blue- 
wing image compared to the in-line image. This is not, how¬ 
ever, clearly seen for some of the events. There exist pixels 
whose co-spatial locations may visually show strong signals 
in both V“ ing and Vn ne images, while the signals are actually 
much larger in the V“ ing than in the Vn ne . The yellow arrows 
on the upper panels mark one of the detected events of the lat¬ 
ter case. The red arrows indicate two examples of events that 
look like detected MUEs (with V“ ing > Viine) but are not be¬ 
cause they do not meet the |V“ ing | > 7cr no i se criterion. We note 
again that in addition to the signal threshold of 7cr no i se for the 
lowing I’ we on ly consider events whose sizes are larger than 
10 pixels (i.e., the same order of the resolution limit of the 
observations). Therefore, all relatively small MUEs, which 
are mostly found in internetwork areas, are excluded in this 
paper. These relatively small events may not be clearly seen 
by visual inspections in Figure 2. 

Similarly, we also identify MDEs from far-red-wing Stokes 
V images and investigate correlations between their spatial lo¬ 
cations with those of MUEs detected from the same scan. The 
bottom-left panel of Figure 2 displays both MUEs and MDEs, 


marked with blue and red circles, respectively. The MDEs are 
derived from the far red wing of the same Stokes V profile 
from which the MUEs are derived. No clear correlation be¬ 
tween the spatial locations of MUEs and MDEs is observed. 

A visual inspection of the spatial locations of both MUEs 
and MDEs on Stokes I continuum images revealed that almost 
all MUEs are located on the granules, close to their edges, 
whereas MDEs are preferably observed on intergranular ar¬ 
eas. This is, however, not so clear for a few of each of the two 
kinds of events whose relatively large areas overlap with the 
edges of granules (i.e., apparently are located on the bound¬ 
aries between granules and intergranules). The bottom right 
panel of Figure 2 represents the Stokes I continuum image 
corresponding to the Stokes V images. Positions of the de¬ 
tected MUEs and MDEs are depicted by small blue and red 
filled circles, respectively (for the sake of simplicity, sizes of 
the circles do not represent the actual sizes of the events). 

A total of 45 MUEs and 47 MDEs are found (using our 
difference approach) in the example scan represented in Fig¬ 
ure 2. We note that our detection algorithm prevents false 
detection (due to, e.g., spatially unresolved and/or relatively 
weak Stokes V signals affected by noise) at the expense of a 
smaller number of detected MUEs/MDEs. 

A much smaller number of MUEs (i.e., 17) are detected 
in the same scan by employing the quotient method (when a 
2noise threshold of |V“ ing | is applied prior to the detection). 

4. ANALYSIS AND STATISTICS 

We employ the same procedure as described in Section 3.3 
(i.e., the same detection algorithm and criteria) to all of the 
image sequences of the four observed passbands. The de¬ 
tected MUEs are then tracked in the time series of images 
and their temporal properties, namely horizontal velocity and 
lifetime, are determined. We use the same code as described 
in Jafarzadeh et al. (2013) for tracking the MUEs. Sizes and 
amplitudes of Stokes V signals of all detected MUEs are also 
calculated. For the latter, we take maximum Stokes V signals 
of all individual pixels within the detected MUEs. We also 
investigate how the MUEs are distributed over network and 
internetwork areas. In addition, we classify Stokes V profiles 
of the detected MUEs. The MUEs are also detected from the 
quotient method (Section 3.2), but only for comparison rea¬ 
sons (i.e., to compare the number densities of the events from 
the two approaches). We do not study MDEs in detail since 
that is beyond the scope of this paper. Their detections are 
only used to investigate possible correlations between their 
spatial locations and those of MUEs. 

To avoid the effect of particular limitations of a particular 
data set on the final quantities (resulting in probable misinter¬ 
pretations), the above analyses are carried out on the data sets 
whose temporal and/or spectral samplings are appropriate for 
certain parameters. Thus, the dynamical properties (lifetime 
and horizontal velocity) of the MUEs are only determined for 
our relatively long time series of images recorded in the 6173, 
6301, and 6302 A lines. The classification of Stokes V profiles 
are only done for the lines whose far blue wings are highly 
(spectrally) sampled, i.e., the 5250 and 6301/2 A lines. The 
other physical parameters (Stokes V signal, size, and line-of- 
sight velocity) are measured for all the datasets. We note that 
the physical and dynamical parameters of the MUEs observed 
in the 6301 and 6302 A pair lines are determined indepen¬ 
dently for the individual lines (with different characteristics, 
e.g., heights of formation; see Table 1). However, simulta- 
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Figure 2. : Examples of far-blue-wing Stokes V (top left), in-line Stokes V (top right), far-red-wing Stokes V (bottom left) and 
Stokes I continuum (bottom right) images from the Fe I 6302 A line. The Stokes V images are saturated to the same scale. 
The blue and red circles mark detected magnetic up flow/downflow events (MUEs/MDEs), respectively. Both events are of either 
polarities and absolute values of their signals are only used for detection purposes. On the upper panels: one of the MUEs is 
magnified, the red arrows indicate two examples of apparent MUE candidates that do not meet signal threshold criterion, and the 
yellow arrows mark an example where the MUE visually shows strong signals in both images. 
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neously observed MUEs are counted only once in the num¬ 
ber statistics. The region (network/intemetwork) association 
of the MUEs is primarily investigated for the latter data set 
(6301 and 6302 A pair lines), from which simultaneous ob¬ 
servations in the Ca II H passband (used to identify network 
areas) are also available. This is, however, also inspected for 
all the lines in our data sets, but from time-averaged Stokes V 
images. 

4.1. Physical and Dynamical Properties 

We detected a total number of 26,517 individual MUEs 
(covering 722,798 pixels in total) from all individual frames 
of the image sequences recorded in the four photospheric 
Fe I passbands (see Section 2 for a description of the data 
sets). It comprises 5150 independent MUEs (when each MUE 
is counted once during its entire lifetime). We note that 
MUEs detected in 6301/2 A lines were counted once since 
the same events were observed in the pair lines. We found 
the same order of magnitude but different rates of occurrence 
of the MUEs detected in the different passbands: 1.8 x 10 -2 , 
3.9 x 10 -2 , and 1.5 x 10 -2 arcsec -2 at any given time from 
images sampled in 5250, 6173, and 6301/2 A lines, respec¬ 
tively. It obviously makes the number of detected events in 
the 6173 A images about twice as larger as those found in the 
other two passbands. 

The area of individual MUEs (i.e., the number of pixels 
that each event contains) and the maximum Stokes V signal 
of each individual pixel within the MUEs were determined. 
We also measured the LOS velocity at each pixel (within the 
MUEs) from bisector positions. The latter are mid-points of 
horizontal segments between the two wings of Stokes I pro¬ 
files. These were measured at different line depths between 
intensity levels of 0% (line core) and 80% (close to contin¬ 
uum level). The difference between the average of the bisec¬ 
tor positions over the different line depths of a profile and a 
reference bisector position results in a mean line shift that is 
translated to the Doppler velocity of that profile (pixel). The 
reference bisector position (or the reference wavelength) was 
computed by averaging mean bisector positions of all profiles 
over the entire FOV. 

We further tracked the MUEs in the image sequences and 
measured their lifetimes and horizontal velocities. The lat¬ 
ter were obtained as a result of dividing their frame-to-frame 
displacements by the time differences between the two con¬ 
secutive frames. 

Distributions of the Stokes V signals (normalized to Stokes 
I continuum), sizes, bisector LOS velocities, lifetimes, and 
horizontal velocities of the MUEs are represented in Fig¬ 
ure 3(a)-(e), respectively. The green, red, yellow, and blue 
lines represent the histograms from 5250, 6173, 6301, and 
6302 A lines, respectively, where applicable. 

Table 2 summarizes ranges and mean values of all param¬ 
eters determined from the different passbands along with the 
number of detected MUEs in each data set. We note that the 
lower and/or upper limits of the parameters presented in Ta¬ 
ble 2 often lie outside the range of the histograms shown in 
Figure 3. All quantities have somewhat different but over¬ 
lapping histograms obtained from observations in different 
wavelengths. Interestingly, the mean values of their Stokes 
V signals, sizes, LOS velocities, and lifetimes (averaged over 
all MUEs) are almost identical to different passbands. The 
upflows are, however, slightly larger for the MUEs observed 
in 6302 A than those in the other three lines. About 1% of the 







Figure 3. : Distributions of Stokes V signals, normalized to 
Stokes I continuum (a), lifetimes (b), line-of-sight (LOS) ve¬ 
locities (c), sizes (d), and horizontal velocities (e) of the de¬ 
tected magnetic upflow events (MUEs). The LOS velocities 
were constructed from bisectors at the wings of the Stokes I 
profiles. The green, red, yellow, and blue histograms repre¬ 
sent the parameters of MUEs observed in the Fe I 5250 A, 
Fe I 6173 A, Fe I 6301 A, and Fe I 6302 A passbands, respec¬ 
tively, where applicable (see the main text). The histograms 
are normalized to their maximum densities. 
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Table 2: Summary of Properties of Magnetic Upflow Events (MUEs) Detected in Different Time Series of Images Sampled in 
Different Photospheric Passbands 




5250 A 

6173 A 

6301 A 

6302 A 

Parameter 

Quantity 

Stokes V/k [%] 

Range 

0.7 , 3.8 

0.6,4.9 

0.6 , 3.6 

0.7,4.4 


Mean 

1.3 

1.3 

1.1 

1.3 

Size [Mm 2 ] 

Range 

0.02,0.25 

0.02,0.45 

0.02,0.36 

0.02,0.40 


Mean 

0.06 

0.05 

0.05 

0.06 

LOS velocity a [km s -1 ] 

Range 

-5.8, 1.2 

-3.0,1.4 

-3.1,2.5 

-3.3, 1.2 


Mean 

-0.8 

-0.8 

-0.7 

-1.1 

Lifetime [s] 

Range 

- 

14,980 

28,1120 

28,1120 


Mean 

- 

113 

114 

113 

Horizontal velocity [km s -1 ] 

Range 

- 

0.0,6.3 

0.0,7.8 

0.0 , 3.2 


Mean 

- 

2.6 

2.1 

1.5 

Number of MUEs 

individual 

2958 

19268 

4291 

4291 


independent 

406 

3730 

1014 

1014 

Rate of occurrence [arcsec -2 ] 


1.8 x 10 -2 

3.9 x icr 2 

1.5 x icr 2 

1.5 x icr 2 


a From bisectors at the wings of Stokes I profiles. 

5250 A MUEs tend to rather have rapid upflows with speeds 
larger than 5 km s -1 . 

There are some differences between the mean values (and 
also upper ranges) of horizontal velocities of the MUEs. In 
particular, the MUEs observed in the 6173 A passband move, 
on average, faster in the horizontal plane compared to those 
found in the 6301/2 A images. While the 6301 MUEs have 
also a lager proper motion, on average, than those observed 
in the 6302 images, the horizontal velocity of a number of 
MUEs detected in the 6173 A and 6301 A datasets reach to a 
speed comparable to that of sound in the solar photosphere. 

The lower limits of the sizes and lifetimes are limited to our 
size criterion of 10 pixels and the cadence of observations, 
respectively. For the latter case, the lifetime of a MUE was 
assumed equal to the cadence of the corresponding dataset 
when it was observed in only one frame. 

4.2. Region association 

The image sequences of 6301/2 A pair lines include simul¬ 
taneous observations of Ca II H filtergrams. Taking into ac¬ 
count the fact that regions with enhanced brightness in Ca II H 
images represent network patches, we create a mask from 
their time average from which the network and internetwork 
areas can also be distinguished in the (spatially and tempo¬ 
rally) co-aligned 6301/2 A images. Using this mask, we in¬ 
vestigate the distribution of MUEs over the network and in¬ 
ternetwork regions. A similar comparison is feasible with 
the temporal average of Stokes V maps where concentrations 
of the magnetic fields representing the network areas are ob¬ 
served. Figure 4 includes all detected MUEs from the entire 
time-series of images observed in 6301/2 A passband, over¬ 
plotted in one frame. The yellow filled contour indicates the 
network areas. Hence, the white regions represent the inter¬ 
network. It turned out that about 58% of all selected MUEs 
are located in internetwork regions, while the rest (42%) are 
observed in network areas. We note that although these frac¬ 
tions are biased by the manual creation of the mask, it re¬ 
veals that the MUEs are approximately equally distributed 
in the network and internetwork areas. Similar results were 
obtained for the MUEs observed in all three lines where net- 
work/intemetwork masks were created from their correspond¬ 
ing (saturated) mean Stokes V images (averaged over their 
entire image sequences). 

4.3. Classification of Stokes V profiles 



0 10 20 30 40 50 

[arcsec] 


Figure 4. : All detected magnetic upflow events (MUEs) over 
the 6301/2 A image sequence of 47 minutes. The yellow and 
white areas represent network and internetwork regions, re¬ 
spectively, determined from the temporal average of simulta¬ 
neously recorded Ca II H images. 

The Stokes V profiles in the MUEs (from 5250 and 
6301/2 A datasets) represent various “anomalous” shapes. 
We carry out a classification of these profiles similar to a 
strategy described by Sigwarth (2001). The V profiles in the 
6173 A MUEs were poorly sampled in the far blue wing (re¬ 
sulting in unclear shapes), thus they are excluded from this 
classification. 

Regardless of any assumptions for their origins and of the 
fact that the different shapes of Stokes V profiles may cor¬ 
respond to different types of events, we classify them solely 
based on their deviations from a normal (double-lobed, anti¬ 
symmetric) Stokes V profile. This is conducted by visual in¬ 
spections of the mean profiles (averaged over all pixels of the 
individual events) according to their appearances. The ratio 
between amplitudes of blue and red lobes are used to subdi¬ 
vide asymmetric double-lobed Stokes V profiles (see below). 
Only “lobes” lying above the corresponding 3 a noise levels 
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are considered as real signals, hence are included in the clas¬ 
sification. 

We find four classes of Stokes V profiles as follows (of 
which one includes two sub-classes): 

(I) Asymmetric double-lobed profiles have two asymmet¬ 
ric lobes with opposite polarities. Their zero-crossing 
points are shifted toward the blue wavelengths in re¬ 
spect to the reference position (i.e., the rest wavelength 
of line core of corresponding lines; see Table 1). 

(a) Profiles with small asymmetries have a ratio be¬ 
tween their blue and red lobes that is smaller than 
50%. 

(b) Profiles with large asymmetries have one lobe 
with an amplitude larger than the other by a factor 
of at least two. 

(II) Asymmetric single-lobed profiles have only one lobe 
whose amplitude lies above the 3 a noise level. 

(III) Double-humped profiles have two lobes of the same po¬ 
larity. 

(IV) Double-lobed profiles with an extra blue-shifted bump 
consist an extra, highly blue shifted bump that is of op¬ 
posite polarity to the main double-lobed profile. These 
profiles look similar to those of typical Stokes Q and U, 
shifted toward blue wavelengths. 

For the sake of consistency, we used same terminology as 
in the literature (e.g., Grossmann-Doerth et al. 2000; Sigwarth 
2001; Borrero et al. 2013) where possible. 

Figure 5 illustrates examples of the different classes of the 
observed profiles along with a typical example of a normal 
(reference) profile (bottom right). The latter represents the 
V profile of a network pixel. In each sub-figure, three im¬ 
ages are also plotted: the far-blue-wing Stokes V (top), the 
in-line Stokes V (middle), and Stokes I continuum (bottom). 
In all three images, the green dotted-dashed circles mark the 
detected MUE. For comparisons, all the MUEs (including the 
one included in the green circles) and MDEs are outlined by 
blue and red contours on the continuum image. The Stokes V 
profiles of all pixels within the detected MUE are plotted in 
gray in each panel. The mean V profiles are overlaid in black 
along with their sampled wavelength positions (blue filled cir¬ 
cles). The vertical lines mark the rest wavelengths of line core 
of the 6301 and 6302 A spectra. 

Among all V profiles of the 6301 /2 A pair lines, the fraction 
of different classes is (I) 27%, (II) 42%, (III) 23%, and (IV) 
3%. The remaining profiles (~ 5%) were excluded from our 
classification because of their peculiar or unclear shapes (i.e., 
they showed more complex appearances). These fractions are 
26%, 39%, 17%, and 5% for classes (I) to (IV), respectively, 
from the V profiles of the 5250 A line. In addition, about 1% 
of the 5250 A Stokes V profiles appeared to be nearly normal 
(perfectly antisymmetric). Twelve percent of the profiles of 
the latter line falls into the unclear cases. 

5. DISCUSSION AND CONCLUSIONS 

We have presented a thorough observational study of MUEs 
from full Stokes observations of four photospheric magnet¬ 
ically sensitive lines recorded at high spatial, temporal, and 
spectral resolutions with the S ST/CRISP instrument. We have 


proposed a robust definition from which the MUEs have been 
detected. The physical and dynamical properties of the MUEs 
have been provided and their Stokes V profiles have been clas¬ 
sified. 

Our detection method (based on the difference between ab¬ 
solute values of significant Stokes V signals in far-blue-wing 
and those in in-line wavelength positions; Section 3.1), to¬ 
gether with the high quality and largely sampled profiles, re¬ 
vealed a larger number of MUEs at any given time (larger 
by one to two orders of magnitude; 2.0 x 10 -2 arcsec -2 
MUEs covering « 0.25% of the entire quiet-Sun area, on 
average) compared to those reported in earlier studies (« 
2.0 x 10 -3 arc sec -2 in Borrero et al. 2010 and Martinez Pillet 
et al. 2011; « 7.0-9.0 x 10 -4 arcsec -2 in Borrero et al. 2013 
and Rubio da Costa et al. 2015). Some of the previous works 
were, however, limited by, e.g., few wavelength positions in 
their observed spectra and/or they were particularly aimed for 
enormous blueshifts in Stokes V signals. These factors could 
result in a relatively small number of detected events. In the 
present work, by contrast, we aimed to study all Stokes V 
signals which were blueshifted from the inner flanks of the 
studied lines. We note that our number of detected events was 
even influenced by our restrictive identification criteria with 
which we excluded MUEs smaller than 10 pixels in size and 
those with Stokes V signals smaller than 7cr no i se . These have, 
however, reduced the chance of false detections. 

Borrero et al. (2010) found their events by identifying sig¬ 
nificant Stokes V signals (signals larger than « 4cr no i se ), in 
a single continuum position at -227 mA from the line core 
of Fe I 5250.653 A. Any significant Stokes V signal was in¬ 
terpreted to stem from strongly blueshifted Fe I 5250.653 A 
components. The corresponding Doppler speed of 12 kms -1 
(to the wavelength difference of 227 mA) lead Borrero et al. 
(2010) to claim discovery of supersonic MUEs. 

Presented in Figure 6 is an example where we find signif¬ 
icant signals (signals larger than 4cr no i se « 0.64% I c ) in all 
five wavelength positions in the far blue wing of 6302 A line 
(panels (a)-(e)). Panel (f) is the averaged in-line positions, 
given as the reference (where no signal is presented at the 
spatial location of the MUE). The five wavelength positions 
(the same as those marked by blue circles with minus signs 
in Figure 1(d)) were sampled in -307, -263, -219, -175, and 
-131 mA from the line core. These correspond to Doppler ve¬ 
locities of-14.6, -12.5, -10.4, -8.3, and -6.2 kms _1 , respec¬ 
tively. Thus, one could conclude that any of these rapid (su¬ 
personic) velocities is associated with the blueshifted Stokes 
V signals, depending on the availability of these wavelength 
positions in the employed data. This may consequently re¬ 
sult in a doubtful interpretation of Stokes V signals at a single 
wavelength positions in the far blue wing in terms of super¬ 
sonic upflow events. We note that this event is associated with 
a bisector LOS velocity of -2.1 kms -1 . We also note that all 
the detected MUEs do not show signals in all single wave¬ 
length positions of the far blue wings, but in some cases sig¬ 
nals only exist in the bluest points whereas in other cases they 
exist only in the reddest positions. In any case, there is signif¬ 
icant signal in, at least one of the five wavelength positions (or 
similarly for the other lines). Hence, the integrated far wings 
also provides the advantage of detecting a larger number of 
events compared to an approach based on one wavelength po¬ 
sition in the far wings (for event detection). 

Although Borrero et al. (2013) also formed in-line and far- 
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Figure 5. : Examples of different shapes (classes) of Stokes V profiles of the magnetic upflow events (MUEs): (la) asymmetric 
but clear double-lobed profile, (lb) asymmetric with a weak second lobe, (II) largely asymmetric, single-lobed profile, (III) a 
profile with two same-polarity lobes, and (IV) double-lobed profile with an extra blueshifted bump. The bottom right panel 
represents a typical example of a normal (double-lobe, antisymmetric) profile of a non-MUE as a reference. Plotted in gray 
are the profiles corresponding to all individual pixels of the MUEs marked by the green circles on the right panels. The mean 
profiles (black lines) as well as the sampled wavelength positions (blue filled circles) are overlaid. The vertical lines mark the 
rest wavelengths of line cores of the 6301 and 6302A spectra. The red dash-dotted lines indicate the 3cr noise level of the 
corresponding averaged profiles. The right panels to the profiles show the images of far-blue-wing Stokes V (top), in-line Stokes 
V (middle), and continuum Stokes I (bottom), with side lengths of ^ 2 arcsec. The in-line V images are saturated to the same 
scale as their corresponding far-blue-wing V images. The blue and red contours in the continuum images indicate the locations 
of MUEs and MDEs in the field-of-view, respectively. The yellow cross in the bottom-right panel mark a network pixel for which 
the reference profile is plotted. 


blue-wing images from their coarsely sampled profiles, their 
lowing I actually represented very far blue-wing (or continuum- 
only) positions. Hence, that lead them to find the smallest 
number of events among the other studies (122 events, cov¬ 
ering 857 pixels in a total area of 1.3 x 10 5 arcsec 2 ). How¬ 


ever, we note that their study was aimed at only searching for 
highly blueshifted V signals. Here, we formed images with a 
relatively high S/N by averaging several wavelength positions 
at the wings and about the cores of the sampled lines prior to 
event detection. Hence, even with a more conservative signal 
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Figure 6.: One example of a magnetic upflow event (MUE) representing significant signals in all five far-blue-wing wavelength 
positions of a 6302 A Stokes V profile (panels (a)-(e)). The numbers in the upper left corners of the panels (a)-(e) show the 
sampling wavelength positions in mA from the line core (the blue circles with minus signs in Figure 1(d)), which correspond 
to Doppler velocities of -14.6, -12.5, -10.4, -8.3, and -6.2 kms -1 , respectively. The blue contours outline the event where 
S/N > 4. The corresponding in-line Stokes V image is shown in panel (f). All color bars show the Stokes V signals in percentage, 
normalized to the Stokes I continuum. 


threshold of 7cr no i se of the treated images, we found a larger 
number of events also using the quotient method. 

In order to investigate the influence of detection methods on 
the number of identified events, we also employed the quo¬ 
tient approach of Borrero et al. (2013) (quotient of absolute 
values of the far-blue-wing and in-line signals; Section 3.2) on 
the same data sets used in the present work. Implementing the 
quotient method using the same criteria as in our difference 
approach (i.e., |V“ ing | > 7cr no i se and size> 10 pixels) results in 
a smaller number of MUEs, by a factor of « 2.5 on average, 
compared to those we found earlier. With the same criterion 
as in Borrero et al. (2013) (i.e., when |V“ ing /Viine| > 4), we 
find thousands of pixels meeting this single condition, where 
no signal threshold is considered prior to or after the detec¬ 
tion. The reason for this is that in the absence of a signal 
threshold, pixels with very weak signals (comparable with the 
noise level) could also satisfy |V“ ing /Viine| > 4. 

The largest MUEs cover an area of about 0.45 Mm 2 (« 
1 arcsec 2 ), while they are on average 0.06 Mm 2 in size. 
These are larger by a factor of two on average, compared 
to those found by Borrero et al. (2010) and Rubio da Costa 
et al. (2015). The Stokes V signals of the MUEs range be¬ 
tween « 0.6-4.0 x 10 -2 I c for all the passbands (reaching to 
a maximum signal of 4.9 x 10 -2 I c in some extreme cases) 
with a mean value of 1.3 x 10 -2 I c . The MUEs live, on aver¬ 
age, for « 110 s, however, their lifetime distributions are ex¬ 
tended toward 10 minutes (and a little bit beyond for a small 
fraction of the events). This is closely in agreement with 
the mean values of 81 s and 86 s reported by Borrero et al. 
(2010) and Rubio da Costa et al. (2015), respectively. The 
LOS velocities of our detected MUEs have a mean value of 
-0.7 to -1.1 kms -1 that is consistent with those from Rubio 
da Costa et al. (2015). The average of sizes, Stokes V sig¬ 
nals, LOS velocities, and lifetimes of the MUEs are almost 
identical in observations in different wavelengths (the life¬ 
times were not determined for 5250 A MUEs because of the 
relatively short image-sequences of those observations). The 
6173 A MUEs appeared to have a larger proper motion (mean 
value of 2.6 kms -1 ) compared to those detected in the 6301 
and 6302 A images (moving with a mean speed of 2.1, and 
1.5 kms -1 , respectively). Those from 6173 and 6301 A tend 
to have relatively wide distributions of horizontal velocity ex¬ 


tending toward the sound speed in the photosphere. This is not 
surprising, though, since these events are found to be mostly 
located at the edge of granules where rapid (sometimes super¬ 
sonic) horizontal convective flows exist (e.g., Cattaneo et al. 
1989; Solanki et al. 1996; Rybak et al. 2004; Bellot Rubio 
2009; Nordlund et al. 2009; Vitas et al. 2011). 

A larger number of MUEs (by a factor of two) were also 
identified in 6173 A far-blue-wing images compared to those 
found in the other passbands. 

The differences between the properties of the MUEs ob¬ 
served in the different passbands as well as the wide distribu¬ 
tions of their physical quantities could depend on the intrin¬ 
sic strength of the lines and/or on their excitation energies. 
These result in different heights of formation and in different 
temperature sensitivities of the lines. Thus, the variation of 
the physical parameters of the MUEs could be dependent on 
heights where they are placed. These differences could also 
be due to, e.g., different conditions of the different data sets 
observed in different times and at different regions (such as 
data quality, S/N, and level of solar activity). 

MUEs smaller than 10 pixels in area and those with V sig¬ 
nals lower than 7cr no j se were excluded in this study. These 
events may have different physical and dynamical properties. 
Hence, their actual distributions could be biased by these se¬ 
lection criteria which have, however, secured our results from 
the influence of falsely detected unresolved and/or relatively 
weak events. 

We found about the same fraction of MUEs in network and 
internetwork areas. The spatial locations of the MUEs did 
not show any clear correlation with those of MDEs. There 
are, however, individual cases where an MUE is apparently 
located in the immediate vicinity of an MDE. 

Furthermore, we classified Stokes V profiles associated 
with the detected MUEs based on their appearances. We 
found them to fall into four classes: asymmetric double-lobed 
profiles (class I), asymmetric single-lobed profiles (class II), 
double-humped profiles (class III), and double-lobed pro¬ 
files with an extra bump (class IV). The latter had an extra 
blueshifted bump in addition to a double-lobed profile (similar 
to those in class I), visually comparable with the appearance 
of Stokes Q or U profiles. 

Different scenarios have been proposed for the physi¬ 
cal processes involved in producing MUEs associated with 
single-lobed V profiles. By inverting a larger spectrally sam- 
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pled line of IMaX Fe I 5250 A (compared to that of Bor- 
rero et al. 2010), Borrero et al. (2013) found coexistence 
of blue- and redshifted flows in about half of their detected 
events. From that, and also from observations of inclined 
magnetic fields of opposite polarities in the immediate vicin¬ 
ity of the events, they concluded that magnetic reconnection 
could be a possible physical mechanism for producing the 
extremely blueshifted magnetized flows. The latter process 
was confirmed by Quintero Noda et al. (2013) who showed 
that interactions between emerging granular-scale loops and 
pre-existing fields in intergranular regions of SUNRlSE/IMaX 
data could produce the strong upflows detected earlier. Later, 
Quintero Noda et al. (2014) found that their detected MUEs 
from Hinode/SP images were often accompanied by MDEs, 
inferring structures corresponding to D-loops. They con¬ 
cluded that siphon flows along arched, magnetic flux tubes 
could explain those field configurations (including the mag¬ 
netic flows along them). 

MHD simulations of emergence and cancellation of mag¬ 
netic flux in the solar photosphere led Danilovic et al. (2015) 
to conclude that magnetic flux emergence is enough to pro¬ 
duce the rapid magnetic flows where magnetic reconnection 
may occur but not necessarily. 

Note that all above mechanisms discuss the magnetic up¬ 
flows corresponding to single-lobed Stokes V profiles. Differ¬ 
ent physical processes may be involved in producing MUEs 
with different shapes of Stokes V profiles (as we found in the 
present work). 

The highly asymmetric profiles (e.g., the single-lobed pro¬ 
files; the class II) are, in general, interpreted as the results of 
strong gradients in the LOS velocity and/or the vector mag¬ 
netic field (e.g., Illing et al. 1975; Landi DeglTnnocenti & 
Landolfi 1983; Grossmann-Doerth et al. 2000; see also Sainz 
Dalda et al. 2012). This may also include the profiles in 
classes la and lb whose second lobes have stronger signals 
compared to those in class II. The level of observed asymme¬ 
try in V profiles, as a result of velocity gradients, has been 
shown to be dependent on the strength (magnetic sensitiv¬ 
ity) of the line, i.e., the velocity gradient produces almost no 
asymmetry in weak lines (Solanki & Pahlke 1988). This may 
also, however, depend on the S/N of the observations. In ad¬ 
dition to the line strength, the dependence of the asymmetries 
have been described on Zeeman splitting, Doppler shift, and 
line width (Grossmann-Doerth et al. 1989). The dependency 
of the line asymmetry on the line strength and the Doppler 
shift could possibly explain the somewhat different fractions 
of asymmetric V profiles we found from the 5250 and the 
6301/2 Alines. 

The double-humped profiles, class (III), as well as the 
double-lobed profiles with an extra blueshifted bump, class 
(IV), could be examples of mixed polarity Stokes V pro¬ 
files. In this case, these profiles could be produced as a re¬ 
sult of, e.g., a spatially unresolved mixture of opposite polar¬ 
ities (Grossmann-Doerth et al. 2000; Sigwarth 2001). Socas- 
Navarro & Manso Sainz (2005) found a few supersonic MUEs 
with three-lobed Stokes V profiles, similar to those that fell 
into our class (IV). They described the extra blueshifted bump 
as a possible signature of an exploding magnetic element due 
to aborted convective collapse (Bellot Rubio et al. 2001). The 
two classes of Stokes V profiles (i.e., the double-humped and 
the double-lobed with an extra bump) could, however, be pro¬ 
duced from a region of unipolar magnetic field (using a two- 
layer model atmosphere) by Steiner (2000), who also pro¬ 


vided a wide variety of Stokes V profiles. They showed that 
while one lobe of the double-humped profiles (as a particular 
case of g-like Stokes V profiles) was in absorption, the other 
lobe could be in partial or complete emission. The latter was 
explained as a result of temperature inversion where the tem¬ 
perature of the magnetic layer surpassed the temperature of 
the line-core-forming regions. Thus, Steiner (2000) showed 
that the line depth plays an important role in forming these 
types of V profiles. 

We note that large blueshifted signals in Stokes V pro¬ 
files do not imply the presence of large (supersonic) veloci¬ 
ties from our bisector analysis; even so, they cannot be com¬ 
pletely ruled out either. In general, the largely Doppler-shifted 
signals could be due to, e.g., high temperatures that could 
increase the Doppler width, a large line shift as a result of 
strong flows, a very high pressure causing strong broaden¬ 
ing of the line, or a turbulence velocity making an additional 
broadening. Therefore, at least for a part of the flows, the 
large blueshifted signals could correspond to very large ve¬ 
locities. Determining these large velocities may, however, not 
be straightforward. The latter needs careful interpretations of 
the profiles, e.g., from Stokes inversions and/or from MHD 
simulations to distinguish the possible thermal widths from 
the fast flows. 

To summarize, we found a larger number of MUEs that was 
larger by one to two orders of magnitude than previously re¬ 
ported. This has been shown to be dependent on detection ap¬ 
proach, on S/N, and possibly on wavelength and resolutions 
of the employed data sets. The MUEs happened to be approx¬ 
imately equally distributed in network and internetwork areas. 
A larger number of events (by a factor of « 2) were identified 
in 6173 A data sets compared to those observed in 5250 A 
and 6301/2 A. The MUEs detected in the former images tend 
to move faster (by a factor of 1.7 on average) in the horizontal 
plane than those found in 6302 A spectra. The MUEs repre¬ 
sent different shapes of Stokes V profiles of which only less 
than half of them are single-lobed. The rest of the profiles ap¬ 
pear to be either asymmetric double-lobed or illustrate signa¬ 
tures of mixed polarities. Some debates also exist on the ori¬ 
gins and on the natures of MUEs (of different kinds). Hence, 
their tentative interpretations need further investigations from, 
e.g., detailed Stokes inversions and/or realistic MHD simula¬ 
tions of various magnetically sensitive lines. It also includes 
a study of their structures and temporal evolutions, character¬ 
ization of the different classes of V profiles associated with 
MUEs, and chasing them in the higher atmospheric layers of 
the Sun. 
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